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0^ , its radius of convergence leading to a pre-asymptotic exponential decay. The 

time of crossover from an exponential to a power law is also predicted. The 
theoretical predictions are confirmed by numerical simulation. Applications 
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0$ , to conductance fluctuations in quantum dots are discussed. 
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I. INTRODUCTION 

Exponential distribution of escape times from chaotic scattering systems should be ex- 
pected only if the associated repeller is hyperbolic. For intermittent repellers one expects 
asymptotic power law decay fl],[|. Nevertheless, in numerical simulations one often observes 
what appears to be a perfect exponential and an expected crossover to a power law 

may be hard to detect, because it may occur after a long time where it is difficult to obtain 
descent statistics. 
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The importance of intermittency cannot be overemphasized. A generic Hamiltonian 
system exhibits a mixed phase space structure. A typical trajectory is intermittently trapped 
close to the stable islands ||. But even fully chaotic billiards may exhibit intermittency, 
typically if they have neutrally stable orbits. Popular billiards such as the Stadium and the 
Sinai billiards are of this type. 

A quantum dot is an open scattering system in two dimensions, obtained by connecting 
leads to a cavity. Inspired by Quantum Chaos research, one likes to contrast shapes of 
the cavity corresponding to chaotic motion, like the stadium, with shapes corresponding 
to integrable motion, like the rectangle or the square. Both extreme cases are sensitive to 
naturally occurring imperfections and one naturally ends up with mixed phase space systems 
where one component hopefully dominates. Consequently, the signals of underlying chaos 
or integrability don't show up as clear cut as one would have hoped for. 

Much of the analysis of these problems has been numerical and heuristic. To strengthen 
the theoretical analysis we will, in this paper, apply periodic orbit theory and cycle ex- 
pansions to make quantitative predictions concerning the asymptotic distribution of escape 
times from an intermittent map. We will demonstrate that a pre- asymptotic exponential 
escape law is associated with a pair of complex conjugate zeroes of the zeta function beyond 
its domain of convergence. This zero will be computed with a simple resummation technique 
0. The truly asymptotic escape distribution will be a power law, and is associated with a 
branch point of the zeta function. The strength of this power law will also be provided by 
the resummation scheme, whereas the particular power is known from analytic argument for 
our particular model system. The relative magnitudes of the pre-exponential and the power 
law will yield a good estimate of the cross over time, which will be surprisingly high. 

II. ESCAPE AND PERIODIC ORBITS 

Much of the early work on cycle expansion |7| was concerned with escape from (hyper- 
bolic) repellers, so we can follow ref. |F] rather closely when deriving the basic formulas 
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relating escape to the periodic orbits of the repeller. 

Consider a 1-d map, on some interval /, with L monotone branches fi(x) where 1 < i < L. 
Each branch fi(x) is defined on an interval Jj. A generating partition is then given by 
£« = {hj 2 ... I L }. We want the map to admit an unrestricted symbolic dynamics. We 
therefore require all branches to map their domain fi(Ii) = I onto some interval I D C^ 1 ' 
covering C^ l >. A trajectory is considered to escape when some iterate of the map x 

The n'th level partition = {I g ; n q = n} can be constructed iteratively. Here q are 
words of length n built from the alphabet A = {i; < i < L}. An interval is thus defined 
recursively according to 

h q = fr\Q , (i) 

where iq is the concatenation of letter i with word q. A concrete example will be given in 
eq. (|T6|) and fig. [I]. Next define the characteristic function for the n'th level partition 

x w W = E^W . ( 2 ) 

q 

where 

!1 X G I a 
(3) 
X £ Ig 

An initial point surviving n iterations must be contained in C^ n \ Starting from an initial 
(normalized) distribution po(x) we can express the fraction that survives n iterations as 

r„= / Po (x)x {n \x)dx . (4) 



We choose the distribution po(x) to be uniform on the interval /. The survival probability 
is then given by 

(n) 

T n = aJ2\ I q\ , ( 5 ) 

q 

where 
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a" 1 = / dx = \I\ . (6) 



1 

Assuming hyperbolicity the size of I q can be related to the stability A q = -£;f n (%) \xeq 
of periodic orbit q according to 

where b q = 0(\I\), can be bounded close to the size of I. This results from the fact that 
f n (I q ) = I, the smallness of \I q \ and the fact that derivatives can be bounded due to 
hyperbolicity. We will eventually relax the assumption of hyperbolicity, but for the moment 
we'll stick to it. 

The survival fraction can now be bounded by the periodic orbit sum according to 

(n) , (n) 

C iW T,TT1<^<C 2 (N) (8) 

q l^ql q \ ly q\ 

for all n > N . For large N (and assuming hyperbolicity) C\(N) and C2(N) can be chosen 
close to unity. 

The periodic orbit sum in @ will be denoted Z n 

(n) * 

and can be rewritten as a sum over primitive periodic orbits (period n p ) and their repetitions 

Zn = T,n P t 6 iff ■ (10) 

p r =l l A J?l 

It is closely related to the trace of the Perron- Frobenius operator 

r 00 5 

J p r =l x l 
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FIG. 1. The intermittent map (|l6| ) for the parameter values s = 0.7 and p = 1.2. The map is 
defined on the interval I. Below the map is also shown the partitions = {1$, 1%}, and . 
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By introducing the zeta function 



^ni 1 -^) • (12) 



the periodic orbit sum Z n can be expressed as a contour integral 



Z n = zr— / Z 



]o S C l (z)jdz , (13) 



2ni Jy \dz 

where the small contour 7 encircles the origin in negative direction. 

The expansion of the zeta function to a power series is usually referred to as a cycle 
expansion. 

C\z)=Y,C n Z n . (14) 

This representation converges up to the leading singularity, unlike the product representation 
(O) which diverges at (nontrivial) zeroes. If the zeta function (~ l {z) is analytic in a disk 
extending beyond the leading zero zq, then the periodic orbit sum Z n , and hence the survival 
probability T n , will decay asymptotically as 

Z n ~ V = e-" n (15) 

where k = log z is the escape rate. 

We will introduce intermittency in connection with a specific model. We then consider 
an intermittent map x t— > f(x) with two branches L = 2, and where / is chosen as the unit 
interval: 

{fo(x) = x (1 +p (2x) s ) x E I = \x; < x < £ } 
" , (16) 

/i(x) = 2a: - 1 xeh = {x]l/2<x<\} 

The map is intermittent if s > and allows escape if p > 1. The map is shown in fig. 
[I], together with some of its partitions. The right edge of the left branch I , denoted £, is 
implicitly defined by /o(£) = 1- The particle is considered to escape if £(s) < x < |. 

The intermittent property is related to the fact that the cycle is neutrally stable 
/'(0) = 1. Consequently, cycle stabilities can no longer be exponentially bounded with 
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length. This loss of hyperbolicity makes it difficult to relate the survival probability T n to 
the periodic orbit sums Z n . Indeed, eq. (0) is brutally violated in some cases, as can been 
realized from the following example. The problem of intermittency is best represented by 
the family of periodic orbits 10 fc . It follows from eqs. ([1]) and (|T6| ) that |iiQfc| = ||/o fc l- It 
can be shown H that I If) k I ~ TTT7 and thus 



l io* 



k i/ s 

which should be compared with the asymptotic behavior of the stabilities 



(17) 



VAio* ~ ^y s (18) 



1.25- 




The difference in power laws seems to spoil every possibility of a bound like (|8|). However, 
eq. (0) is not necessary for that purpose. It suffices if the ratio 

IA I np 

R P = — H \W (19) 
n p k=i 

stays bounded. That it to say that it suffices if the average size of the intervals along a cycle 
can be related to the stability, rather than each interval separately. Here S denotes the shift 
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operator. We check this numerically on two sequences: 10 fc and 110 fc , the former being most 
sensitive to intermittency. The result is plotted in fig. [| We note that for both sequences, 
R p appear to tend to well defined limits, where 10 fc lead to the largest deviation from unity. 
Indeed, it is reasonable to assume that the sequence 10 fc provides a lower bound 

R p > lim R 1Qk Vp . (20) 

k— >oo 

In view of this, the numerical results strongly suggest that R p stays bounded, and that, for 
this particular system, C\ in eq. (H) can be chosen as Ci = 0.5 and C2 presumably close to 
unity. This is a surprisingly low price to pay for the complication of intermittency. 

The sizes of the intervals Ion has no relation whatsoever to the stability of the cycle 0, 
which is unity. We exclude the intervals Jo n from our considerations by pruning the fixpoint 
from the zeta function 

^-SO-S) ■ (21) 

The contribution from I n to T n can be added separately if required. 

Since the result rely on summation along periodic orbits, it might break down for some 
choices of the initial distribution po(x) were such a summation is not carried out. 



III. RESUMMATION AND SIMULATION 

After having argued that the survival probability T n still can be bounded close to periodic 
orbit sums z n we turn to the problem of computing the asymptotics of these periodic orbit 
sums. The coefficients of the cycle expansion flTID for the map flT6|) decay asymptotically as 



1 (22) 



^ n l+l/8 ' 

which induces a singularity of the type {l-z) 1 ' 8 in the zeta function |J . If 1/s is an integer, 
the singularity is (1 — z) 1 ^ log(l — z). 

To evaluate the periodic orbit sum it is convenient to consider a resummation of the zeta 
function around the branch point z — 1. 
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C\z)=^c i z i = ^a i {l-z) i + {l-z) x l s< >rb i {l-z) i . (23) 

i=0 i=Q i=Q 

In practical calculations one has only a finite number of coefficients q, < % < n c of the 
cycle expansion at disposal. Here n c is the cutoff in (topological) length. In || we proposed 
a simple resummation scheme for the computation of the coefficients cij and 6« in (p3|). We 
replace the infinite in (^) sums by finite sums of increasing degrees, n a and rib, and require 
that 

n a n b n c 

^a i (l-z) i + (l-z) 1 / s £b i (l-z) i = '£c n z n + 0(z n+1 ) . (24) 

i=0 i=0 n=0 

One expands (1 - z) i( - +1 /^ in binomial sums (series). If n a + rib + 2 = n + 1 this leads to a 
solvable linear system of equations yielding the coefficients <2j and 6j. It is natural to require 
that \n a + - — rib\ < 1 so that the maximal powers of the two sums in eq. (p4|) are adjacent. 

If the zeta function is entire in the entire z plane (except for the branch cut) the periodic 
orbit sum can be written 



(25) 



The sum is over all zeroes z a of the zeta function (assuming they are not degenerated) and 
the contour 7 CMt goes round the branch cut in positive direction. If poles and/or natural 
boundaries are present, expression ( p5[ ) must be accordingly modified. 

The leading asymptotic behavior is provided by the vicinity of the branch point z = 1, 
and is found to be || 

M l l _ 

a s 1 (1 — l/s) n 1 / 5 

The relevant ratio &o/ a o> obtained from the resummation scheme, versus cutoff length n c 
is plotted in fig. |3|. In all numerical work we have used the parameters s = 0.7 and p = 1.2. 
and computed all periodic orbits up to length 20. 
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FIG. 3. The ratio bo/ao versus cutoff length n c 



There is also a pair of complex conjugate zeroes, z = x ± iy close to the branch cut. 
They contribute both to the sum over zeroes and to the integral around the cut in (p5|) . 
But since their imaginary part yo is small, they will, in effect, contribute a factor Xq n to the 
periodic orbit sum Z n . 

This zero will dominate Z n in some range C n < n cross before the asymptotic power 
law sets in. 

In fig. |] we study the convergence of the real and imaginary part of zq obtained from 
the resummation scheme above, for different cutoffs n c . The zero is computed by Newton- 



Raphson iteration of the left hand side of (|24|), after the resummation has been done. Again 
we note that the analytic continuation technique works quite satisfactorily. 
The probability of escaping at iteration n is 



Pn 



n-l 



(27) 



We get for this distribution 



Pn 



Xq 
b 1 



a s 2 r(l-l/s) nW/s lb cros 



10 



< n < n cross 



(28) 



Here we have neglected the interval Jo™ > having the same asymptotic decay law as the periodic 
orbit sum Z n . Due to the uncertainty in the bounds (§) it can be neglected. 

The crossover n = n cross takes place when the two terms in fl2"8|) are of comparable 
magnitude. For our standard set of parameters {p = 1.2, s = 0.7) it is found to be n cross « 
300. 




FIG. 4. The real and imaginary part of log zq versus cutoff length n c . 

The check our predictions we run a simulation of the system. The result can be seen in 
fig. ([|). We note that the slope of the exponential, the power and magnitude of the power 
law, as well as the crossover time agrees very well with our predictions. 

A reader still in any doubt on the effectiveness of cycle expansions should consider the 
following. The simulation in fig. |5] averaged over 10 8 initial point, yet, in itself the result 
would not very conclusive. A direct evaluation of the periodic orbit sum up to say n = 600 
would require roughly 10 170 periodic orbits. We have not bothered to perform such a cross 
check! But a resummed cycle expansion provides reliable answers with a length cutoff as 
low as n = 15, corresponding to 4719 prime cycles. Admittedly, we benefitted from knowing 
the asymptotic power law of the cycle expansion. However, if this is not the situation, this 
power law is easily extracted if one uses stability ordering ||. 
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n 

FIG. 5. Distribution of escape times obtained from simulation (shaky curve) and the pre expo- 
nential (full line) and the asymptotic power law (dashed line) obtained from resummation. 



IV. A ZERO IS CUT IN TWO 

The occurrence of a dominating zero beyond the branch point is, in fact, very natural. 
Consider the one-parameter family of zeta functions. 

= 5(1-^) • (») 

For j3 small enough there is a leading zero z (/3) within the domain of convergence \zo(/3) | < 1. 
This is related to the topological pressure PJlQf according to V{(3) = — logZo(P)- F° r 



instance V(0) is the topological entropy. For a certain f3 (actually the fractal dimension of 
the repeller) the zero collides with the branch point z — 1, splits into two, and continue to 
move out beyond the branch point. 

In fig. ^ we plot the logarithm of the leading zero (— log z (j3)) versus (3. It is obtained 
from a resummation analogous to the one discussed above, cf. ||. It can be interpreted as 
the topological pressure only as long it is real. 
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FIG. 6. The real and imaginary part of the leading zero versus the thermodynamic parameter (3 



V. MESOSCOPIC DISCUSSION 



The particular form of the distribution of escape times does depend on the initial distri- 
bution po(x). In this paper we have restricted ourselves to a uniform initial distribution. To 
model chaotic scattering one must imagine that particles can be injected according to any 
distribution. For example, one can construct a chaotic scatterer from a bounded billiard by 
drilling holes wherever on the boundary and injecting particles from different angles. This 



may even effect the asymptotic power law ||11|| . Periodic orbit theories can also account for 
other initial distributions than uniform. However, the preceding discussions about relat- 
ing periodic orbit sums to survival probabilities warns us to be cautious when doing so for 
intermittent systems. 

As it appears, the general rule of thumb, first an exponential, then a cross over to some 
powerlaw, can be extended to open Hamiltonian systems with a mixed phase space structure 

ill- 

An immediate application concerns conductance fluctuations in quantum dots |j!2| . The 
Fourier transform C(x) of the correlation function C(Ak) = {T(k)T(k + Ak))k, where T{k) 



13 



is the transmission as function of the Fermi wave number, can, after several approximations, 
be related to the escape distribution p(L) ||T2f 



poo 

C(x) ~ / dLp(L + x)p{L) . (30) 

J 

If there is a crossover to a power law in p(L) there will be an associated crossover in C(x). 
For an intermittent chaotic systems, the crossover time may be very long - the quasi regular 
region component of phase space will not make itself noticed until very long times. If the 
elastic mean free path of electrons is much shorter than the length corresponding to the 
crossover time, the quasi regular component will never ve detected in this type of exper- 
iments. Or the other way around, a small deviation from an integrable structure induces 
chaotic layers in phase space. This chaotic layers may lead to exponential escape for small 
times, and the experimental outcome may very well resemble predictions for fully chaotic 
systems. 

In experiments a (weak) magnetic field is a more natural control parameter than the 
Fermi energy. Instead of the distribution of dwelling times one has to consider the distribu- 
tion of enclosed area, a related but more subtle concept which we plan to address in future 
work. One has observed Lorentzians shape (predicted for chaotic systems ) of the so called 
weak localization peak even in near integrable structure [[DJ. This has been attributed to 



naturally occurring imperfections [fn,|l5l and rhymes well with the classical considerations 
above. Admittedly, we have now moved far from our original intermittent map and entered 
the realm of speculations. What we do want to point out in this letter is that these kind of 
problem are well suited for periodic orbit computations - zeta functions is a powerful tool for 
making long time predictions, even for intermittent chaos, once the problems of analytical 
continuation can be overcome. 
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